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2 ; Abstract 
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Q ■ The phase spaces of the two- and three-frequency sine-Gordon models are exam- 

■^ ! ined in the framework of truncated conformal space approach. The focus is mainly 

on a tricritical point in the phase space of the three-frequency model. We give sub- 
stantial evidence that this point exists. We also find the critical line in the phase 
_ , space and present TCSA data showing the change of the spectrum on the critical 

Q I line as the tricritical endpoint is approached. We find a few points of the line of first 

i-C ' order transition as well. 
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1 Introduction 

The sine-Gordon model in (1 -|- 1) dimensions has attracted interest long time ago for 
the reason that it appears in several areas of physics, nevertheless it is an integrable field 
theory, that can be used to study non-perturbative quantum field theory phenomena. The 
areas of application include statistical mechanics of one-dimensional quantum spin chains 
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and nonlinear optics among many others — see the introduction of [^ for a representative 
hst with references. 

In the present paper we investigate a non-integrable extension of the sine-Gordon 
model called multi-frequency sine-Gordon model, in which the scalar potential consists 
of several cosine terms with different frequencies. It is suggested in fT| that this model 
can be used to give more refined approximation to some of the physical situations where 
the ordinary sine-Gordon model can be used. A feature of the multi-frequency model 
that is new compared to the usual sine-Gordon model is — apart form non-integrability 
— that phase transition can occur as the coupling constants are tuned. We concentrate 
our attention to this property. It should be noted that such a phase transition is related 
to the evolution of the particle spectrum of the theory as the coupling constants vary, 
and we shall use the massgap and other characteristics of the energy spectrum as order 
parameter. 

Our investigation is a continuation of the work done on the double-frequency case in 
[3]. We use the truncated conformal space approach (TCSA), which is a non-perturbative 
numerical method to compute the lowest levels of the spectrum of the Hamiltonian opera- 
tor of perturbed conformal field theories in finite volume. The applicability and reliability 
of this method was thoroughly investigated in j3|, and it was shown that the existence, 
nature and location of the phase transition can be established by this method, although 
rather large number of dimensions are needed for satisfactory precision. In particular, the 
existence and location of an Ising type transition was established in the double-frequency 
model (DSG) for the ratio 1/2 of the frequencies, verifying a prediction by [l| based on 
perturbation theory and classical arguments. We extend these investigations to the ratio 
1/3 and to the three-frequency model (at the ratio 1/2/3 of the frequencies), in which a 
tricritical point and first order transition are expected to be found. The numerical nature 
of the TCSA makes it necessary to choose specific values for the frequencies. 

After introducing the multi-frequency sine-Gordon model and describing basic proper- 
ties of it in Section O we briefly review the TCSA framework for the model (which is the 
framework described in jH] adapted to the multi-frequency model) in Section 01 In Section 
m we give a description of the phase structure of the classical two- and three-frequency 
model, which serves as a reference for the investigations in quantum theory. The n- 
frequency case is also considered briefly. Exact and elementary analytic methods can be 
applied to the classical case, and the results are more general than in the quantum case. 
Section ini is devoted to theoretical considerations on the signatures of 1st and 2nd order 
phase transitions in the framework of perturbed conformal held theory in flnite volume. 
Most of these considerations, which are necessary for the evaluation of the TCSA data, 
can also be found in 0. Finally, in Section Inland [3 we present the results we obtained by 
TCSA on the phase structure of the two- and three-frequency model, which are the main 
results in the paper. A summary of the results and conclusions are given in Section |H1 



2 The multi- frequency sine- Gordon model 

2.1 The definition and basic properties of the model 

The action of the multi-frequency sine-Gordon model (MSG) is 

Amsg = jdtjdx Qa^^a'^^ - V {<!>)] , 
where 

n 

V^($) = ^/i,cos(A$ + <5,) 

i 

is the potential, which contains n cosine terms. $ is a real scalar field defined on the 
two-dimensional Minkowski space R^+^, /3j G M are the frequencies, /3j 7^ Pj Hi ^ j, fii 
are the coupling constants (of dimension mass^ at the classical level) and 5j G M are the 
phases in the terms of the potential. 

Two cases can be distinguished according to the periodicity properties of the potential. 
The first one is the rational case, when V"($) is a trigonometric polynomial: the ratios 
of the frequencies Pi are rational and the potential is periodic. Let the period of the 
potential be 27ir in this case. The target space of the field $ can be compactified: 

$ = $ + 2kr7r, 

where A; G N can be chosen arbitrarily. The model obtained in this way is called the 
/c-folded MSG. The well-known classical sine-Gordon model corresponds to n = 1, A; = 1. 
The other case is the irrational one, when the potential is not periodic and no such 
folding can be made. The irrational case is much more complicated than the rational 
one, so we restrict our attention to the rational case in the present paper. We remark 
here only that although V{^) always has a finite infimum, it does not necessarily admit 
an absolute minimum. The potential V{^) can always be written uniquely as a sum 
V{^) = Vi{^) + V^2($) + ••• + Vfc($), where the terms Vi,...,Vk are periodic but any 
sum of any of these terms is not periodic. V"($) has an absolute minimum if and only 
if Vi(<l>), ..., Vfc($) have a common absolute minimum. This occurs for special choice of 
the Si, if the values of Pi are given. In particular, if Pi/Pj are irrational for all i j^ j and 
Hi < for all z, then V(^) has a absolute minimum if and only if -f — -5^ = ^^ ^ is 

Pi Pj Pi Pj 

satisfied with some numbers bi G Z, which is equivalent to the case Si = for all i. See 
and |Tj and IH] for further remarks on the irrational case. 

At the quantum level the theory can be regarded as a perturbed conformal field theory: 

Amsg = Ac ft + Apert, 
where 



AcFT= I dt I dx-a^$a^$, 



which is the action of the free scalar particle of zero mass, and 

Apert = fdtfdx{-Vm = -l fdt fdx ^(/i.e^'^'Vft + /i,e-*'^'Vlft), 

where V^ denotes the vertex operator 

V =■ e*'^* • 
which is a primary field with conformal dimensions 

in 






in the unperturbed (conformal) field theory. The upper index ± corresponds to the 
left/right conformal algebra and : : denotes the conformal normal ordering. The dimen- 
sions of the couplings at the quantum level are 

[/i,] = (mass)2-2^% A, = A^^. 

The perturbing operators are relevant only if 

Pi < Svr, (1) 

we restrict ourselves to this case. We also assume that 

A' < 4vr, 

which is a necessary and sufficient condition for the model to be free from ultraviolet 
divergencies in the perturbed conformal field theory framework |HIII1IS1- 

The model has a massgap in general, and it is clear that phase transitions occur in the 
classical version of the model as the coupling constants are tuned (assuming that n> 1). 
It is also expected that there are topologically charged solutions/states in the model |T|. 
In the present paper we shall investigate the sector with zero topological charge, which 
is sufficient for our purposes. We also restrict ourselves to 1-folded models (/c = 1), as it 
is natural to expect that in infinite volume a folding number k ^ 1 results simply in a 
/c-fold multiplication of the spectrum corresponding to k = 1. 

3 The Truncated Conformal Space Approach 

The matrix elements of the total Hamiltonian operator H = Hcft + Hp^rt can be cal- 
culated explicitly (and exactly) between any two elements of the Hilbert space TCcft of 
the conformal field theory. Restricting to a finite dimensional subspace Ho of Hcft by 



introducing an upper conformal energy cutoff H becomes a finite matrix, which can be 
diagonalized numerically to get an approximation of the spectrum. As this description 
shows, TCSA is a variational method. 

The /c-folded free boson as a conformal field theory contains the following primary 
fields: 

Vp^p{z,z) =: explipcpcFriz) + ip(t)cFT{z)] : 

2 -2 

with conformal dimensions A+ = |-, A~ = f-, where p = - + 2Tcrm, p = - — 2'Krm, 
n,m E Z, 

^cft{x, t) = (pcFxix -t) + (J)cft{x + t). 

TicFT is spanned by the states |p, p) = lim^^^^o ^,^(-2? ^)|0) (|0,0) = |0)) and a„j...ami...|p,p), 
where a„^ and 0,^^ are creating operators of Fourier modes on the "conformal plane". The 
conformal generators Lq and Lq are diagonal in this basis. The basis of Tio is obtained by 
taking those elements \v) of the basis above which satisfy the truncation condition 

{v\^Hcft\v) 



{v\\v) 

Ccut is the dimensionless upper conformal energy cutoff and L is the volume of space. We 
restrict ourselves to the sector with zero topological charge {p = p), this being the sector 
containing the ground state (s) and the relevant information for the problem treated in 
this work. In the sector containing the ground state(s) the Lorentz spin of all the states 
is also zero, so we also impose the condition (Lq — Lq)\v) = 0. (The operator Lq — Lq 
commutes with H and He ft as well) . We remark that Ccut serves as an ultraviolet cutoff. 
The matrix elements of if between two elements \a) and \b) of the basis oiTicFT above 
are given by 

H\ 2tt f ^ - c\ 



jafcOAa-Aa.Afc-Aj, 



9-7r / A/f \ ^-' 1^3 

+ — Yl ''9n{fXj)nj l-j^j _^ e-'^^{V-(s„~(3,{l, l))a6<5A„-A„,A,-A„ 

where M is a mass scale of the theory given below, / = LM is the dimensionless volume, 
Xj = 2 — 2Aj, Aj = Ap.; Aq, Aa, Afo, A;, are the conformal weights of the states \a) and 
1 6), c is the central charge of the conformal theory (c = 1 in the present case), and we 
have made a replacement corresponding to 



The "interpolating" mass scale M is 



J' 



where 

are the dimensionless coupling constants (of which only n — 1 are independent). © implies 
that rjj G [0, 1], J2jVj = 1- M depends smoothly on the rjj-s. The precise expression for 
K, is not essential for our problem, we need only that n depends on A only and that it is 
dimensionless. Following [H] we used the formula of f?j: 




The formula Q is written in terms of dimensionless quantities, and the volume (/) 
dependence of H/M is also explicit. It is clear from Q that TCSA gives an exact result 
if / — i> (assuming {TJ) and this limit of the theory is the conformal theory (the massless 
free boson), and the accuracy of the TCSA spectrum decreases at fixed Ccut as / — ;> oo. 
For very large values of / the /-dependence of the spectrum of the TCSA Hamiltonian 
is power-like, as it is determined directly by the l^^~^ coefficients in Q. The TCSA 
Hamiltonian cannot be considered as good approximation for these values oil. Increasing 
the value of Ccut generally increases accuracy, and the TCSA is expected to give better 
approximation for the lower lying energy levels than for the higher ones. 

We denote the (dimensionless) energy levels oi H/M in volume / by ej(/), z = 0, 1, 2, ..., 
and Co < ei < 62 < ... if not stated otherwise. We shall draw conclusions about the 
spectrum at / = oo from the behaviour of the functions €«(/) for low values of i and 
moderately large values of /. 

We refer to |3l IH E] for more details on the TCSA framework. 

4 Phase structure in the classical limit 

4.1 Phase structure of the two- frequency model in the classical 
limit 

The Lagrangian density takes the following form in the two-frequency case: 



1 
2' 



C = -d^<^d^<^ - fi cos(/3<l>) - A cos(a<l> + 6), 



where 

a m 
n and m are coprimes (and the folding number equals to one). 

Proposition 1: Assume that /z, A 7^ 0. Then the following three cases can be distin- 
guished: 

a.) If the function V{^) is symmetric with respect to the reflection $ 1-^ 2$o — ^, 
where $0 is a suitable constant, and n,m > 1, then V has two absolute minima, which 
are mapped into each other by the reflection. We remark that it depends on the value of 
6 whether V has this symmetry or not, and it is easy to give a criterion for the existence 
of this symmetry in terms of n, m and 6. 

b.) If V is symmetric with respect to a reflection as in case a., but n = 1 or m = 1, 
then V has one or two absolute minima depending on the values of /i and A. In this case, 
assuming that n = 1, V can be brought to the form 

V{(^) = -|/i| cos(/?$) + |A| cos(m/3$) 

by an appropriate shift of $. V has two absolute minima if | A//i| > 1/m^, and one absolute 
minimum if |A//i| < 1/m'^. The two absolute minima are mapped into each other by the 
reflection. The second derivative of V is nonzero at the minima if |A//x| 7^ 1/m^, but it 
is zero if |A//i| = 1/m^. In the latter case, the fourth derivative of V at the minimum is 
nonzero. The two minima of V merge and the value of the second derivatives of V at the 
two minima tends to zero as |A//i| approaches l/rri^ from above. 

c.) If V does not satisfy the requirements of a) and b), then V has a single absolute 
minimum. 

We omit the proof of this proposition, which is elementary, although long and not 
completely trivial because of the arbitrariness of n and m. 

If /i or A equals to zero, then V is periodic and has m or n absolute minima, respec- 
tively. See [4J for a detailed investigation of these (integrable) limiting cases. 

The phases of the classical model are determined by the behaviour of absolute minima 
of V{^) as the value of the coupling constants vary. In particular. Proposition 1 implies 
that the phase structure of the two-frequency model is the following: 

The model exhibits an Ising-type second order phase transition at the critical value 

m 

Vc 



1 + m 



of the dimensionless coupling constant rj = \/\jA/i^/\jA + \/\M) ii n = 1 and V has the 
Z2-symmetry introduced above. This critical point separates two massive phases with 
unbroken and spontaneously broken Z2-symmetry. Equivalent statement can be made if 
m = 1. If K is not symmetric, then there is only one massive phase with nondegenerate 
ground state. If m, n 7^ 1 and V is symmetric, then there is one massive phase with 
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doubly degenerate ground state (i.e. the reflection symmetry is spontaneously broken). 
In the limiting cases r] = and ?] = 1 the model is massive and has spontaneously (and 
completely) broken Z„ or Z^ symmetry. 

4.2 Phase structure of the three-frequency model in the classical 
limit 

A complete description of the behaviour of the absolute minima of V for all values of the 
parameters becomes excessively difficult in the three- and higher-frequency cases, so we 
restrict our attention to particular values. The potential in the three-frequency case is 

V{<^) = /ii COs(/3i$) + /i2 cos(/32$ + ^2) + /is cos(/33<^ + Ss). (4) 

We choose the frequency ratios 3:2:1, i.e. 

Pi =P, P2 = Ip, P3 = IP- 

This three-frequency model has a tricritical point if and only if ^2 = ^3 = (and also in a 
few equivalent cases), in this case V is symmetric with respect to the reflection $ \-^ — $. 
In the tricritical point the absolute minimum of V can be located only at or vr. The two 
cases are equivalent, we consider the case when the location of the absolute minimum is 
0. The tricritical point in this case is located at 

^ - _i /^ - ^ 

/i2 6' /is 15 

In this point V^^^O) ^ 0. (The upper index ^^^ denotes the sixth derivative with respect 
to $.) V"*^^^(0) > requires /ii,/i3 < and /i2 > 0. We restrict ourselves to this domain 
and to the values ^2 = ^s = 0. 

The phase diagram is shown in Figure [H The points of the diagram correspond to the 
values of the pair {rji,rj2) of dimensionless parameters. The allowed values constitute the 
left lower triangle, the straight line joining (0, 1) and (1, 0) corresponds to 773 = 0. The 
tricritical point is denoted by t, it is located at 

1 Vq 



1 + Ve + v^' 1 + Ve + v^ 



^ (0.1365,0.3345). 

Ki t V has one single and absolute minimum (at $ = 0). Phase transition occurs when 
the lines 5 and 3 shown in the phase diagram are crossed. Second order Ising-type phase 
transition occurs on 5 and flrst order phase transition occurs on 3. The domain A U 
B U F corresponds to a massive Z2-symmetric phase (with unique ground state). The 
domain E U C corresponds to a massive phase with spontaneously broken Z2-symmetry. 
Characteristic shapes of the potential in the various domains and on the various lines of 
the phase diagram can be seen in Figure El Data applying to the quantum case are also 
shown in Figure [H they will be explained in subsequent sections. 
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Figure 1: Phase diagram of the classical three-frequency sine-Gordon model at (i\l (iij (iz = 
3/2/1, ^1 = (^2 = (^3 = 0, /ii, /i3 < and /i2 > 0. The crosses and the square correspond 
to certain quantum theory values described in Section [3 



vc*) 




C 4,E,7 a,5,t 

Figure 2: Characteristic shapes of the potential 

4.3 n-frequency model in the classical limit 

Let us take the n-frequency model with Pi = ip, i = l..n, and 6i = 0. In this case there 
exist unique values of /ij//ii, i = 2..n so that V{x) has a single global minimum at a; = 
and has no other local minima, and V"{0) = 0, V""{0) = 0, ... y^^'^^O) = also hold. 
The values of Hi/ni are determined by the latter equations. The point corresponding to 
these values of /ij//ii is an n-fold multicritical point in the phase space. The neigbourhood 
of this multicritical point contains m-fold multicritical points for any integer < m < n. 
These statements can be proved using well known properties of analytic functions and 
the fact that V^ is a trigonometric polynomial. We omit the details of the proof. 

5 Signatures of 1st and 2nd order phase transitions in 
finite volume 

The considerations in this section apply to quantum field theory. 

The behaviour of the spectrum is governed by the / -^ limiting conformal field 
theory for small values of /, so e„(/) — eo(/) ~ 1//. Massive phases in infinite volume are 
characterized by the existence of a massgap and the behaviour hm;^oo(en(0~6o(0) = C'„, 
where C„ > are constants. C„ = OifO<n<(i and C„ > if n > d, if the ground 
state has d-fold degeneracy in infinite volume. In a phase with spontaneously broken 
symmetry the spectrum is degenerate in the / -^ cxd limit, in finite volume the degeneracy 
is lifted (at least partially) due to tunneling effects. The resulting energy split between 
the degenerate vacua vanishes exponentially as / -^ oo. 
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In the critical points (in infinite volume) the massgap vanishes and the Hilbert space 
contains a sector that corresponds to the conformal field theory specifying the universality 
class of the critical point. We consider this sector in the following discussion. In finite 
but large volume near the critical point this sector of the theory can be regarded as the 
/ — *> oo limiting conformal theory perturbed by some irrelevant and relevant operators. 
The corresponding TCSA Hamiltonian takes the (generic) form 



27r / , rf \ CiR , \-^ 9-^^ 



2-2A 



ii 



where the i}) are the perturbing fields. This picture gives the following volume dependence 
of energy levels (in the first order of conformal perturbation theory) : 

e^{l) - eo(0 = ^(A,\^ + A7^^J + 5^4/l-2^^ (6) 

where A|^ ^ and ^'J^ ^ are the conformal weights of the state '^ in the / -^ oo limiting 

CFT, A% are constants that also depend on the particular energy eigenstate \E'. The 
presence of irrelevant perturbations (1 — 2A^ < —1) is due to the finiteness of the volume, 
whereas the presence of the relevant perturbations (1 — 2A^ > —1) is caused by the 
deviation of the control parameters from the critical value and by the UV cutoff. 

The location of the critical points can be determined using the criterion of vanishing 
massgap. A more precise method that also allows the determination of the universality 
class of the critical point (i.e. the I — > cxd limiting CFT) is the following: We make 
an assumption that the critical point is in a certain universality class. This assumption 
predicts the set of ip-s, the values of Aj"^^ and A^^^, and the values of the A^-s in ^. 
We take leading terms of the series on the r.h.s. of (0) and determine the value of the 
(A]^^^ + A7/j^)-s and of the A%-s by fitting to the TCSA energy data. The magnitude 

of the v4^-s corresponding to the relevant perturbations measures the deviation from the 
critical point, so if the assumption on the universality class is right, then by tuning the 
coupling constants one should be able to find a (critical) value of the control parameter at 
which these v4^-s are small, the TCSA data are described well by © (terms from higher 
orders of perturbation theory can be included if necessary) in a reasonably large interval of 
the values of/, and the values of the (Aj^^ + A//j^)-s obtained from the TCSA data agree 
with the assumption with good precision. The interval where ^ describes the TCSA data 
well is called the scaling region. This region may be (and in fact is) different for different 
energy levels. We remark that it is also possible to make a theoretical prediction for 
60(0, which allows to extract cir from the TCSA data for eo(/) in principle [H]. However, 
experience (0) shows that the accuracy of the TCSA data is not sufficient to determine 
ciR precisely in this way, so we do not attempt to extract cir directly from the TCSA 
data in this paper. 
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a b 

Figure 3: Typical shape of V{^) in the broken and in the unbroken symmetry phase 

In the classical case a first order phase transition occurs when the absolute minimum 
of the potential becomes a relative minimum and a previously relative minimum becomes 
absolute. In the quantum case this phase transition is characterized by the presence of 
"runaway" energy levels with asymptotic behaviour e{l) ~ cl for large / in the neighbor- 
hood of the transition point, where c is a constant that tends to zero as the transition 
point is approached. The multiplicity of the ground state also changes as the transition 
point is passed if the two phases have different symmetry properties. We remark that 
"runaway" energy levels are present in general whenever a model has unstable vacua. 



6 The phase diagram of the two-frequency model in the 
case -^ = 3,0=3 

We assume that A,/i > 0, and use the parameter 



V 



yp 



^x^ + X^P 



instead of rj in this case to conform with (H]. 

The classical model with a/ [3 = |, 5 = | exhibits an Ising type phase transition at 
77 = 37(1 + 3^). 

Considering the quantum case we proceed along the lines of [HI in this section. The 
numerical nature of the TCSA makes it necessary to choose a finite number of values for 
/3 and I at which calculations are done. One should choose as large values for / as possible, 
the / —^ CX3 limit being of interest. However, the accuracy of TCSA decreases as I grows. 
The accuracy can be improved by taking higher Ccutj but this increases the size of TCSA 
Hamiltonian and the time needed for diagonalization. Experience shows that accuracy 
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decreases for values (3 near to v47r (this is the value where UV divergences appear in 
conformal perturbation theory) although the speed of convergence of the spectrum to 
the / -^ oo asymptotic values increases, and the speed of convergence becomes very low 
for values of /3 near to 0. Taking these properties of the TCSA into consideration and 
following [HI we performed calculations at the values P = Sa/vt/T, iy/n/S, Sy/n/d. We also 
note that the accuracy of the TCSA spectra is severely decreased if V has several (local) 
minima. 

Figure Ola andlHlb show the shape of the classical potential in the phases with broken 
and unbroken Z2-symmetry, respectively. Figures IHa-IHg show TCSA spectra obtained 
at /3 = 4:y/n/3 at various values of fj. The TCSA Hilbert space had dimension 3700, the 
first 12 energy levels are shown in the figures. The highest values of / are chosen so that 
the truncation error still be small (the massgap remain constant). However, the effect of 
truncation is perceptible in Figure |3|b for instance. It can be seen that in the domain 
fj < 0.92 the ground states and the first massive states are doubly degenerate. (They are 
triply degenerate at fj = 0.) "Runaway" energy levels (of constant slope) corresponding 
to the single local minimum of the potential can also be seen (especially clearly in Figure 
I31b). In the domain rj > 0.98 the spectra are massive, but the ground state and the first 
massive state are nondegenerate. In the intermediate domain (especially for fj ~ 0.95) 
the structure of the spectrum changes, no massgap and degeneracy can clearly be seen. 
We obtained similar spectra at /3 = 8-yvr/5 and /3 = Sy^/? as well. As we did not 
see "runaway" energy levels that would have signaled first order phase transition in the 
transitional domain of fj we analyzed the data by looking for a second order Ising type 
phase transition at some critical values fjc{(3). 

The Ising model contains three primary fields: the identity with weights (0,0), the e 
with weights (1/2, 1/2), and a with weights (1/16, 1/16). Since the DSG model exhibits 
the Z2-symmetry for all values of fj, the Z2-odd a and its descendants cannot appear 
as perturbations in the Hamiltonian ^. The only relevant field compatible with the 
Z2-symmetry is e (the contribution of the identity cancels in the relative energy levels). 
The presence of a relevant perturbation e in the Hamiltonian leads to a correction B^^, or 
B^i + C^l to e^(/) — eo(0- The term C^/ is of second order in perturbed CFT. The leading 
irrelevant perturbation (compatible with the Z2-symmetry) is the first descendant of e, 
this gives a correction A^/"^ to e^(/) — eo(0 (in first order). Thus we expect that in a 
large but finite volume range, near fJdP), the volume dependence of the energy levels is 
described well by the formula 

e^{l) - eo(/) = y A + AI-' + B, + ai. (7) 

We fitted this function to the lowest energy levels obtained by TCSA and determined the 
"best" fidP) value by tuning fj in the transition region and looking for whether e2(/) — eo(/) 
continues to decrease along the complete I range (lim/^oo(e2(0 ~eo(0) = only at fjciPj), 
and Bi and Cj are as small as possible. The result is shown in Table [H The fitting was 
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State 


D:. 


A, 


B, 


a 


1 = 1 


0.138 ±0.0005 


-2.0 ±0.02 


-0.0046 ± 0.0001 


2.98 ■ 10"^ ± 9 • 10"^ 


1 = 2 


1.00 ±0.01 


-74 ±2 


0.004 ±0.001 


9- 10-6 ±4 -10-^ 



[3 = 8v^/5, r/ = 0.944 



State 


Di 


A. 


B, 


Q 


i = 1 


0.125 ±0.001 


-2.63 ±0.025 


-0.0002 ±0.0001 


3-10-^±l- 10-6 


i = 2 


1.04 ±0.02 


-152 ±6 


0.0014 ±0.0009 


-1.5- 10-5 ±2 -10-6 



j3 = 4v/^/3, fi = 0.955 



State 


Di 


Ai 


B, 


Q 


i = 1 


0.125 ±0.001 


-4.4 ±0.1 


-0.0009 ±0.0001 


4-10-^±5- 10-^ 


i = 2 


1.00 ±0.02 


-252 ± 9 


0.0002 ±0.0004 


4.3- 10-6 ±6 -10-^ 



P = 8v^/7, f] = 0.961 



Table 1: The results of fitting ^ to the first two excited levels for various values of /3 at 
the estimated critical value of fj 

done in the volume ranges / = 10 - 105, / = 55 - 105; / = 10 - 140, / = 100 - 190; 
/ = 20 — 200, / = 150 — 390. The errors presented come from the fitting process and do 
not contain the truncation errors which are generally much larger. 

The first two energy levels above the ground state correspond to the operators a and e in 
the Ising model. These operators have conformal weights A^*" = 1/16 and A^*^ = 1/2, so 
the exact values for Di and D2 are 



Di = 0.125, 



D, 



1. 



The results of the fits agree quite well with this prediction. 

The TCSA data obtained at the estimated values of fjc using a truncated space with 
dimension 4800, 5300 and 5100 are shown in Figure El These figures show energy levels 
multiplied by //(27r) as functions oil. The constant lines corresponding to the Ising model 
values of D^ are also shown in these figures. 

To summarize, evaluating the TCSA data we found that the phase transition is second 
order and Ising type at the values of /? chosen. This statement is true up to the precision 
of the TCSA data, which is about 10-^ — 10-^ for the values of the volume used. We 
remark that it is not possible to distinguish a second order transition from a (very) weakly 
first order transition by TCSA. 
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a: f] = 



b: fi = 0.3 





c: 1] = 0.7 



d: fi = 0.92 





ISO 2O0 250 300 



e: ?7 = 0.94 



f: 77 = 0.96 




g: ?7 = 0.98 

Figure 4: Change of the spectrum as fj varies from to 1 at /9 = 4:^/n/3. The first 
12 energy levels (including the ground level) relative to the ground level are shown as 
functions of /. 
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[ei(/) - eo(0] ■ V2vr at /3 = 80F/5 and 77 = 0.944 



[ei{l) - eo(/)] ■ 1/271 Sit p = 4v^/3 and r) = 0.955 



400 
1 



[ei(/) - eo(/)] ■ 1/2TC at (3 = 80F/7 and fj = 0.961 
Figure 5: TCSA spectra as functions of / at the estimated critical values of fj 
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Figure 6: The phase diagram of the two-frequency model at a/ [3 = 1/2 and a/P = 1/3 



For values of /3 near we expect second order phase transition, because the model is 
semi-classical in this region and so the correction to the classical potential in the effective 
potential is expected to be small. 

A correction to the classical potential in the effective potential of frequency 2(3/3 is 
possible in principle. A correction with this frequency is — unlike in the case of a/ (3 = 
1/2 — relevant for any values of /3, and it can be verified by elementary calculation that 
it may change the order of the transition outside the semiclassical region if its coefficient 
is sufficiently large (see also [21 [HI). However, the accuracy of TCSA did not allow us 
to perform calculations at values about /3^ > Stt/S and to check the nature of phase 
transition in this domain. 

The phase diagram in the (/3, fj) plane based on the data obtained by TCSA can be 
seen in Figure IHl We took into consideration that fic{0) = 3*^/(1 + 3'^) ~ 0.988 is exactly 
known P = being the classical limit, and at /? = a/Stt the term with frequency /? in the 
potential becomes irrelevant and thus for /? — ^ a/Stt the other term and so the symmetric 
phase is expected to dominate, so lim^^^/g^ f7c(/3) = 0. The figure shows the three values 
of fie obtained from the TCSA data and the two values at (3 = and (3 = VSn. The 
continuous line is obtained by fitting an even polynomial (note that fjc{P) = fjc^—P)) to 
these values and is shown in order to get an idea of the phase transition line. The model 
is in the symmetric phase above the line and in the phase with broken symmetry below 
the line. The data obtained by for the a/f3 = 1/2 case are also shown, the dashed line 
is fitted to these data. 
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7 Phase diagram of the three-frequency model 

We take the same values of the parameters of (JH) as in section 14 .21 namely Pi = P, P2 = 
|/3,/?3 = |/3, 52 = ^3 = 0, yUi,/i3 < 0, yU2 > 0, and investigate the quantum model in this 
case. 

7.1 The tricritical point 

The tricritical Ising model contains 6 primary fields with the following conformal weights: 



and 



3 3 \ / 7 7 



80' 80/ ' V16' 16 ' ■ ^ ' 



The fields corresponding to (jH} are even and those corresponding to are odd with 
respect to the parity, so only the fields corresponding to (jH)) and their descendants can 
contribute to the Hamiltonian ^. Thus the volume dependence of the energy levels near 
the tricritical point should be described well for large / by 

e^{l) - eo(/) = y (A+j,, + A,^,,,) + A^r^'^ + B^f'' + ..., 

where only the leading terms are kept. Searching for the tricritical point we fitted the 
function 

^ A + A,l-'-' + Bf-' (10) 

to the data obtained by TCSA for ej(/) — eo(0 near the estimated location of the tricritical 
point. Best fits are shown in Table El (The errors presented come from the fitting process 
and do not contain the truncation errors which are generally much larger.) The fitting 
was done in the volume ranges / = 50 — 230, / = 110 — 230, the dimension of the truncated 
Hilbert space was 13600. The results of the fitting support the existence of a tricritical 
point located (approximately) at rji = 0.163, 772 = 0.3518. The exact values of Di and D2 
in the tricritical Ising model are 

Di = 0.075, D2 = 0.2. 

The numerical results agree quite well with this prediction. The TCSA spectrum obtained 
at the tricritical point is shown in Figure [3 The values of D^ predicted by the tricritical 
Ising model are also shown in the figure. The dashed lines and + signs are used for odd 
parity states, the continuous lines and x signs are used for even parity states. 
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State 


D^ 


A^ 


B^ 


i = 1 


0.074 ±0.004 


-0.0060 ±0.001 


2.8 ■ 10-5 ± 5 ■ 10-*^ 


i = 2 


0.196 ±0.01 


-0.006 ±0.002 


2.4 ■ 10-5 ± 6 ■ 10-*^ 



p = 8v^/7, Til = 0.163, ri2 = 0.3518 



Table 2: The results of fitting (fTTHl to the first two excited levels in the estimated tricritical 
point 
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Figure 7: [ei{l)—eo{l)]-l/2T!- as functions oil obtained by TCSA at /3 = 8-\/7r/7, rji = 0.163, 
r/2 = 0.3518 
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Vi 


V2 


Vi 


V2 


Vi 


m 


Vi 


V2 


0.01 


0.354 


0.06 


0.355 


0.11 


0.359 


0.16 


0.357 


0.02 


0.354 


0.07 


0.3555 


0.12 


0.36 


0.163 


0.3565 


0.03 


0.354 


0.08 


0.356 


0.13 


0.36 






0.04 


0.354 


0.09 


0.3575 


0.14 


0.3595 






0.05 


0.354 


0.1 


0.3585 


0.15 


0.359 







Table 3: Points of the critical line found by TCSA 

We used a modified version of the TCSA code exploiting the Z2-symmetry of the 
model by taking even and odd basis vectors and taking the Hamiltonian on the even and 
odd subspaces separately, which reduces the total time needed for diagonalization and 
thus allows to take higher Ccut values. 

The TCSA data for the first two excited states fit quite well to the prediction of 
the tricritical Ising model, the energy levels of the next two excited states also show 
correspondence with the prediction. Clear correspondence cannot be seen for higher 
levels. 

7.2 The critical line 

The points of the critical line we found using TCSA are listed in Table [HI The value of rji 
was chosen and fixed in advance, and then 772 was estimated in the same way as described 
in Sectional These points are also marked in Figure [l] by crosses. The dimension of the 
truncated Hilbert space was 10269 in these calculations, which corresponded to Ccut = 17. 
The value of (3 was S^tt/S. Figures |Hla-l show the TCSA spectra (especially the lowest 
lying energy levels) obtained in these points as well as the values of Di corresponding to 
both the critical and the tricritical Ising model. Crosses are used for odd parity states, 
squares are used for even parity states. It can be seen that moving on the critical line 
in the phase space towards the tricritical endpoint the finite volume spectrum changes 
continuously. In the r]i < 0.11 domain the spectra (especially the first two levels) corre- 
spond clearly to phase transitions in the Ising universality class. At rji = 0.11 the first 
excited level already appears to correspond to the prediction of the tricritical Ising model, 
whereas the second excited level still has the behaviour predicted by the Ising model. It 
would be very interesting to know the large volume behaviour (and infinite volume limit) 
of the first excited level, but the precision of TCSA does not allow to determine it. What 
we can see is that there is no sign in the TCSA data that the first excited level follows the 
predictions of the Ising model in the large volume limit. In the domain 0.11 < rji < 0.16 
there is a spectacular rearrangement of the higher energy levels (already observable in the 
rji < 0.11 domain), and at r]i = 0.16 the second excited level also appears to correspond 
to the prediction of the tricritical Ising model. We regard therefore the point rji = 0.16, 
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rj2 = 0.357 to be the tricritical endpoint of the critical line, this point is marked by a 
square in Figure QJ (We did not aspire to determine the value of r/i more precisely for 
this value of /?.) Each figure shows the spectrum in the volume interval / = 0..200, and 
for the large values / ~ 200 the truncation error is always conspicuous. 
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f: rii = 0.11 



g: r/i = 0.12 h: r]i = 0.13 
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Figure 8: [ei{l) — eo{l)] ■ 1/2tt as functions of / obtained by TCSA at /? = Sy^/S and at 
various points (771,^72) lying on the critical line, the predictions of the critical Ising model 
(continuous horizontal lines) for D,, the predictions of the tricritical Ising model (dashed 
horizontal lines) for Di 

Moreover, the value of / where the truncation errors become large gets smaller and smaller 
as the tricritical point is approached, which corresponds to the fact that a tricritical point 
as renormalization group fixed point is more repelling than an Ising type one. Increasing 
rji further the behaviour corresponding to the tricritical point would rapidly disappear 
from the finite volume spectrum. 

7.3 The line of first order transition 

Figures[nia-j show TCSA spectra obtained at rji = 0.6 and at various values of 772 between 
0.17 and 0.30. The energy levels are shown as compared to the lowest level. The dimension 
of the truncated Hilbert space was 6597 in these calculations, which corresponded to 
Ccut = 16. The value of /3 was 8-\/vr/5. Dashed lines are used for odd parity levels and 
continuous lines for even parity levels. 

At ri2 = 0.17 (and also for ri2 < 0.17) the ground state is unique for all values of /, 
whereas doubly degenerate runaway levels can also be seen. At 772 = 0.19, however, the 
ground level is doubly degenerate for small values of the volume but becomes nondegen- 
erate in large volume. 
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c: r/2 = 0.19 



d: r/2 = 0.20 






e: r]2 = 0.21 



f: r/2 = 0.22 




g: V2 = 0.23 



h: r/2 = 0.24 
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i: r/2 = 0.25 



j: V2 = 0.30 



Figure 9: [ei{l) — eo(/)], i = 0..19 as functions of / obtained by TCSA at /3 
rji = 0.6 and at various values of 772 



80F/5, 



i.e. there is a value oil where the doubly degenerate runaway level and the single level cross 
each other. The slope of the runaway levels become smaller and smaller as 772 is increased, 
and the crossing point also moves towards larger and larger values. At 772 = 0.3 the 
ground state is already doubly degenerate for all values of / and nondegenerate runaway 
levels are present. These features indicate that a first order phase transition occurs at an 
intermediate value of 772. The behaviour of the finite volume spectra around the transition 
point seen in the figures implies that a precise determination of the transition point from 
the TCSA data in a direct way would require precise data for large values of l. For this 
reason we did not aspire to find many points of the line of first order phase transition, 
we did TCSA calculations only at r/i = 0.4 and r/i = 0.6, and we roughly estimated the 
location of the first order phase transition at these values. Our estimation based on the 
TCSA data are r]2 = 0.30 and 772 = 0.21. These points are also marked in Figure [l] by 
crosses. The FigureslHla-j show the first 20 levels in the domain / = 0..200. Unfortunately 
the truncation effect is large in most of this domain of Z, however we expect that those 
qualitative features of the spectra which we allude to are correct. 

8 Conclusions 

In this work we investigated selected special cases of the multi-frequency sine-Gordon 
theory, especially the structure of their phase space, continuing the work begun in [3j. 
Concerning the classical limit we found that the only possible phase transition in the 
(rational) two-frequency case is a second order Ising-type transition. The three-frequency 
model has some new qualitative features compared to the two-frequency model: a tricriti- 
cal point which is at the end of a critical line can also be found, and first order transition is 
possible as well. The phase space of the n-frequency model contains n-fold critical points, 
otherwise we do not expect new qualitative features compared to the three-frequency 
model. 
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The numerical (TCSA) calculations in the quantum case yielded the following results: 
we found the same type of phase transitions as in the classical case, in particular we were 
able to determine the second order Ising nature and the location of the phase transition 
in the two-frequency model with good precision. Quantum corrections did not alter the 
nature of the phase transition. The accuracy of the TCSA also allowed us to find the 
tricritical point in the three-frequency case, which we regard the main result of the paper. 
This demanded much more numerical work than the two-frequency model for the following 
reasons: the two-frequency model has only one-dimensional phase space, whereas the 
phase space is two-dimensional in the three-frequency model; and the tricritical point as 
a renormalization group fixed point is more repelling than the Ising- type critical point, 
so we had to take larger truncated space. Furthermore, we found several points of the 
critical line and observed how the TCSA spectrum changes as the tricritical endpoint is 
approached. It would be interesting to investigate this with considerably better precision. 
Although the focus was mainly on the tricritical point, we also investigated the first order 
phase transition in the three-frequency model and we found that the first order nature can 
be established by TCSA, but the location of the transition can be determined much less 
accurately than that of the second-order transition. We expect that multicritical points in 
the universality classes of further elements of the discrete unitary series could be found in 
the higher-frequency models, but increasing accuracy is needed, because these multicritical 
points are more and more repelling. Finally, we remark that the quantum corrections did 
not alter the nature of the phase transitions in any of the cases we investigated, and in the 
classical and quantum theory the location of transition points is almost the same. The 
investigation of the multi-frequency model with irrational frequency ratios is still an open 
problem. We also remark that the particle content of the multi-frequency sine-Gordon 
model could also be investigated by the method of semiclassical quantization [H]. 
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